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A  lized-biiit  band  structure  interpolation  scheme 


applied  to  the  rocksalt  struct nr e  con  pounds  TiC,  TiN  and  TiO 


Sehun  Kin  and  R.  Stanley  Williams 

Department  of  Chemistry  and  Biochemistry,  and 
Solid  State  Sciences  Center 
University  of  California,  Los  Angeles,  CA  90024 


A  mixed-basis  band  structure  interpolation  scheme  has  been  applied  to 
the  rocksalt  compounds  TiC,  TiN  and  TiO  by  fitting  to 
first-principles  calculations.  For  TiC  and  TiN,  semi-empirical  band 
structures  were  obtained  by  adjusting  the  interpolation  parameters  to 
improve  the  agreement  between  the  calculated  bands  and  data  f ran 
previous  Angle-Resolved  Photoelectron  Spectroscopy  and  X-ray 
Photoelectron  Spectroscopy  measurements.  The  use  of  the 
eigenfunctions  to  plot  approximate  valence  charge  densities  to 
illustrate  bonding  effects  is  described. 


N  V 


I.  Introduction 


Mixed-basis  band  structure  interpolation  schemes  (MBBSIS),  which 
employ  a  combination  of  d-like  orbitals  and  plane  waves,  have  proven 
to  be  extremely  valuable  computational  tools  for  examining  the 
electronic  structures  of  the  d-band  metals.  Shiith  and  various 
coworkers*-*  have  developed  a  fee  d-band  interpolation  scheme  over 
the  past  decade  and  applied  it  as  an  aid  in  interpreting 
photoemission  and  UV  reflectance  data.  An  extension  of  Smith's 
scheme  was  ntilixed  by  Kim,  et  al.,^  to  examine  the  band  structure  of 
several  fluorite-atructure  compounds  with  d-states  in  the  valence 
band.  In  this  paper,  the  MBBSIS  is  applied  to  compounds  with  the 
rocksalt  structure,  and  the  use  of  the  eigenfunctions  to  plot 
qualitative  valence  charge  densities  is  illustrated. 


The  rocksalt  structure  can  be  described  as  a  fee  Bravais  lattice 
with  a  basis  consisting  of  an  A  atom  at  (0,0,0)  and  a  B  atom  at  the 
position  (a/2) (1,0,0)  of  the  conventional  cubic  cell.  A  large  nnnber 
of  carbides,  nitrides  and  oxides  of  the  transition  metals  crystallize 
in  the  rocksalt  structure.  All  carbides  and  nitrides  of  the  Group  IV 
and  V  transition  metals  have  common  properties:  high  melting  point, 
hardness,  brittleness  and  metallic  conductivity.  All  of  these 
physical  properties  are  related  to  the  electronic  structure  of  the 


refractory  compounds.  The  importance  of  the  electronic  structure  for 
these  compounds  has  motivated  many  theoretical  and  experimental 
investigations.  General  overviews  of  studies  on  the  fourth  and  fifth 


group  transition  metal  carbides,  nitrides  and  oxides  are  given  in 
review  articles^'^ 

This  investigation  was  aotivated  in  part  by  the  observation  that 
the  empirical  pseudopotential  aethod  (EPM)  calculations  of  the  TiC 

o 

electronic  structure  by  Alvard,  et  al.,  who  employed  a  nonlocal 
pseudopotential  to  treat  the  Ti  d-states.  produced  valence  charge 

densities  that  disagreed  severely  with  those  obtained  from 
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first-principles  calculations  and  experimental  investigations. 

The  EPM  calculations  predicted  that  the  local  p-d  charge  density 

maximum  near  the  carbon  atom  (the  C  2s  orbitals  were  not  included  in 

the  calculation)  was  twice  as  large  as  that  near  the  titanium  atom. 

while  the  first  principles  calculations  showed  relatively  comparable 

total  (ie.  including  C  2s)  charge  densities  around  the  two  atoms. 

The  intent  of  this  study  was  to  see  if  the  1BBSIS  could  be  used  to 

obtain  reliable  band  structures  and  also  qualitatively  correct  charge 

densities  for  complex  systems  such  as  TiC,  TiN  and  TiO.  These 

systems  have  been  studied  extensively,  and  thus  there  exists  a  great 

deal  of  data  to  compare  with  the  MBBSIS  results. 

Since  the  rocksalt  structure  has  the  same  space  group  symmetry 
as  fee  crystals,  the  same  MBBSIS  may  be  used  for  both  lattices.  The 
situation  is  similar  to  that  for  the  previous  study  of  the  fluorite 
structure  compounds.^  In  that  investigation,  a  MBBSIS  similar  to  that 
of  Stoith,  et  al.*-*  was  used,  except  that  more  plane  waves  were 
included  in  the  basis  set.^  In  Sec.  II,  the  spplication  of  the  MBBSIS 


to  the  rocksalt  structure  con  pounds  and  the  calculation  of  charge 
densities  is  described.  In  Sec.  III.  the  energy  bands,  the  total 
density-of-states  (TDOS)  and  the  electron  density  plots  calculated  by 
the  IBBSIS  are  presented  and  discussed.  Finally.  Sec.  IV  contains 
the  conclusions  of  this  study. 


II.  Interpolation  Scheme  and  Charge  Density  Calculations 

Our  method  of  calculation  has  been  described  previously.^  It  is 
nearly  identical  to  the  MBBSIS  of  Sbith,*-*  except  that  the  plane 
wave  block  of  the  Hamiltonian  is  a  39x39  submatrix,  instead  of  a 
16x16.  The  additional  plane  waves  are  required  to  represent  the 
wavefunctions  of  the  non-metal  atom  in  the  compound.  The  total  local 
pseudopotential  of  the  rocksalt  compound  AB  can  be  written  as 
follows : 

Vg(AB>  *=  Vg(A)  +  VG(B)cos(fi*t),  (1) 

M  Ml 

where  the  d-bands  originate  from  the  A  atom  and  (a/2) (1,0,0)  is  a 
vector  connecting  the  transition  metal  atom  (A)  to  the  nometal  atom 
(B)  in  the  primitive  unit  cell.  Por  the  G  vectors  of  the  reciprocal 
lattice  of  the  rocksalt  structure,  the  values  of  cos(G*t)  are  limited 
to  +1  and  determined  only  by  the  magnitude  of  g,  not  its  direction. 
Thus,  the  values  of  Vg  for  AB  are  all  constants,  with 

vg(a)  -  vg(b)  for  IgI  =  VT,  Vn 

A* 

Vg(AB)  -  (2) 

vG(A)  +  vg(B)  for  lul  -  o,  2,  VT,  VTT, 

where  |G I  is  in  units  of  2n/a.  As  long  as  the  d-orbital  contribution 
of  the  B  atoms  to  the  band  structure  is  negligible,  the  interpolation 
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Htailtonimn  of  the  rocksalt  structure  is  identical  to  that  of  the  fee 
atructure. 

The  parameters  of  the  MBBSIS  for  a  particular  system  were 
determined  by  applying  a  nonlinear  least  squares  method7*  to  fit  the 
energy  eigenvalues  to  those  of  first-principles  calculations. 74-7* 
The  parameters  usually  converged  to  values  that  yielded  average 
deviations  (up  to  the  tenth  band  of  the  rocksalt  structure  compounds) 
of  about  0.2  eV  after  5~<  iterations.  Soni-empirical  energy  bands 
were  then  generated  by  adjusting  the  parameters  of  the  MBBSIS  to 
improve  the  agreement  with  available  ARPES  data  for  TiC14,17  and 
TiN.7*  After  the  fitting  parameters  for  a  band  structure  were 
determined,  a  TDOS  for  the  system  was  calculated  by  the  method 
described  previously.* 

The  electron  density  associated  with  each  occupied  state  E 

n,  k 

can  be  defined  as  the  square  of  its  eigenfunction: 


(3> 

where  #^^(5)  and  dj  ^(jc)  are  the  39  Orthogonal ized  Plane  Waves  (OBf) 
and  the  5  d-wave  functions,  respectively,  and  b®  tn(j  c°  are  the 
coefficients  of  the  basis  functions  that  result  fr cm  diagonalizing 
the  MBBSIS  Hamiltonian  at  point  k  in  the  BZ.  The  set  of  OlWfs  is 
given  by 
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(4 ) 


where  the  plane  wave  Ik^  -  o_^^ei^  ~,  0  ia  the  primitive  nnit  cell 
volume,  and  •  (1-  Z<dj  In  the  present  charge  density 

calculations,  since  the  contribution  of  the  orthogonal ization  terms 


in  eq.  (4)  was  negligible,  the  were  approximated  as  being 

simply  plane  waves  normalized  to  a  charge  density  of  two  electrons 


per  primitive  unit  cell.  The  five  d-wave  functions  (d  .a  ,d 

xy’yz  'xz ' 

<*x2_y2  and  djx2_r2)  can  be  written  as 


dj.k<S>  “  *3d<*)Y2j(k)' 


(5) 


where  ^3^(5)  is  the  radial  part  of  the  d-wave  functions  and  the 
*2j(k>  sre  properly  normalized  spherical  harmonics.  In  the  MBBSIS 
calculations  employed,  no  explicit  form  for  the  radial  part  of  the 
d-wave  function  was  used.  An  alternative  form  of  the  interpolation 
scheme,  which  incorporated  Slater  radial  d-wave  functions  and  used 
the  orbital  exponent  as  a  fitting  parameter,  as  in  Ref.  5,  was 
tested.  However,  simply  using  Smith’s  form  for  the  Hamiltonian  terms 
of  the  MBBSIS,  then  substituting  the  single-zeta  radial  wave  function 
for  Ti  3d  taken  from  Clement!  and  Roetti's  tables*’5  into  Eq.  3, 
produced  charge  densities  nearly  identical  with  those  calculated 
using  the  alternative  scheme. 


The  electron  density  calculated  as  a  function  of  r  for  an 
arbitrary  k-point  does  not  necessarily  have  the  symmetry  of  the 
crystal.  Imposing  the  crystal  symmetry  on  the  charge  density  can  be 


accomplished  by  either  summing  ^(£)  over  the  star  of  £  or  applying 
the  appropriate  space  group  operations  on  ^(£> .  In  the 

calculations  reported  here,  the  48  space  gronp  operations  of  0^  were 
applied  to  the  charge  density  calculated  at  each  value  of  £.  This 
procedure  was  computationally  faster  than  evaluating  p^  ^<£)  over  the 

star  of  k,  *nd  yielded  the  same  results.  Total  valence  charge 
densities  were  calculated  by  suing  over  all  the  valence  bands  at  10 
special  k-points*®.  A  large  ntnber  of  simplifications  and 
approximations  have  been  used  to  calculate  these  charge  densities. 
The  primary  value  of  such  calculations,  as  long  as  they  are 
qualitatively  correct,  is  to  be  able  to  visualize  the  electronic 
structure  of  a  relatively  complex  system  in  real  space  (ie.  a  bond 
picture  of  the  system)  without  resorting  to  extremely  expensive 
computations  or  requiring  an  enormous  first-principles  basis  of 
eigenvectors. 

All  the  calculations  were  performed  using  a  VAX  11/780  computer. 
The  Central  Processing  Unit  (CFO)  time  of  the  VAX  required  to  produce 
the  eigenvalues  at  a  single  k-point  was  3.7  sec.,  and  a  value  of  the 
charge  density  at  a  particular  r  for  a  band  state  E  .  required  0.3 


III.  Results  and  Discussion 

Since  relativistic  effects  are  mall  for  the  3-d  transition 
aietals.  the  MBBSIS  calculations  for  the  Ti  compounds  were  performed 
without  including  a  spin-orbit  psrameter.  Most  of  the  energy  band 

calculations  for  the  transition  metals  of  group  IV  and  V  reviewed  in 

Refs.  [6,7]  were  limited  to  no  higher  in  energy  than  the  tenth  energy 
band.  Recently,  Linearized  Augmented  Plane  Tave  (LAPIf) 
calculations^ .15  £or  jiC  and  TiN  generated  band  structures  up  to 
35  eV  above  the  Fermi  energy  (Ep} .  The  MBBSIS  was  used  to  fit  these 

LAPW  energy  bands  up  to  the  twelfth  and  part  of  the  thirteenth  bands 

(18  eV  above  Ep)  using  45~50  LAJPV  energy  eigenvalues  at  4~5  high 
symmetry  points  in  the  BZ.  Since  the  LAFW  N  2s  bands  were  emitted  in 
Ref.  15,  the  N  2s  bands  of  the  MBBSIS  were  initially  fitted  to  those 
of  the  A CT  calculations  by  Neckel  et  al.^®  The  average  deviation  from 
the  LAPW  energy  bands  was  ~0.5  eV  for  TiC  and  TiN.  These  values  are 
~3  times  larger  than  the  average  deviation  for  fits  to  the  TiO  band 
structure,  since  the  A JW  calculations  for  this  compound  did  not 

include  eige  . axues  above  the  tenth  band.  The  MBBSIS  for  TiO  was 

fitted  to  the  early  band  structure  of  Ern  and  Switendicfc^  because 
these  authors  listed  their  eigenvalues  in  a  Table  and  thus  made 
possible  a  more  precise  comparison  with  a  previously  calculated  band 
structure.  However,  the  band  structure  plot  of  Ref.  16  agreed  very 

well  with  that  shown  by  Neckel,  et  al.  in  a  more  recent 


The  most  severe  deviations  of  the  MBBSIS  energy  bands  from  the 
LAPW  eigenvalues  occurred  at  the  second  band  of  ("^  and  the  first  band 

of  symmetry,  and  the  largest  contributions  to  the  average 

deviation  come  from  energies  above  the  tenth  band  for  TiC  and  TiN. 
The  agreement  for  the  ten  bands  lowest  in  energy  was  generally  within 
0.2  eV.  In  the  case  of  TiO,  for  which  energy  bands  were  generated 
only  up  to  the  tenth  band  in  Fig.  3,  the  average  deviation  from  the 
AWf  energy  bands***  was  comparable  to  the  results  of  a  fit  to  an 
improved  LCAO  method  that  used  45  disposable  parameters.  Thus,  the 
MBBSIS  is  able  to  fit  the  rocksal t-com pound  bands  as  well  as  the  LCAO 
method  with  fewer  parameters. 

The  TDOS  for  Ti  compounds  calculated  by  the  MBBSIS  agreed  very 
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well  those  with  calculated  from  LCAO  bands  fitted  to  ATO  results 
(see  Figs.  1-3).  Most  of  the  experimental  studies  to  date  have 
provided  information  about  the  TDOS  of  these  compounds,  but  do  not 
contain  k-resolved  information  about  the  band  structures.  In  recent 
studies,  however,  detailed  band  mapping  was  performed  using  ARPES  to 
obtain  information  about  band  dispersions  and  surface  states  for 
TiC**’*^  and  TiN.*^  These  data  allowed  direct  comparison  between 
calculated  band  structures  and  ARPES  results.  Noticeable 
discrepancies  between  the  LAPW  calculations  and  ARPES  data  can  be 
found.  The  second  band  of  A^  symmetry  and  the  first  band  of  A^ 
symmetry  calculated  by  the  LAPW  technique  lie  closer  to  the  Fermi 
level  than  those  determined  from  ARPES  data  for  TiC  and  TiN.  For 
example,  in  the  case  of  TiN,  the  experimentally  determined  position 
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of  the  ri5  point  is  -3.4(+0.1)  eV.  of  the  Xj.  point  -4.0(+0.2)  eV, 
and  of  the  X^,  point  -6.4(+0.2)  eV15.  compared  with  calculated 
▼alnes  of  -2.2,  -3.3,  and  -5.9  eV respectively.  Also,  the  LAW 
energy  bands  for  TiN^  along  f-X  did  not  show  a  gap  between  the  first 
band  of  A^  and  the  first  band  of  symmetry  along  the  A  direction, 
bnt  the  ARPES  results  indicated  a  substantial  gap1*'  (see  Fig.  4). 

The  MBBSIS  bands  of  TiC  and  TiN,  which  were  initially  fitted  to 

14  15 

the  LAW  energy  bands,  *  were  adjusted  to  agree  with  the 
experimental  results.  '  The  resulting  band  structures  and  TDOS 

plots  are  shown  in  Fig.  1  and  2  for  TiC  and  TiN,  respectively,  and 

the  MBBSIS  results  for  TiO  are  shown  in  Fig.  3.  The  empirically 

adjusted  energy  bands  for  TiN  agree  very  well  with  the  ARPES  data 
along  f-X  and  show  clear  a  band  gap  between  the  third  (A^)  and  fourth 
(A^,)  bands,  as  shown  in  Fig.  4. 

The  possible  caveat  to  adjusting  the  valence  bands  to  agree  with 

photoemission  data  is  that  such  bands  may  not  reflect  the  true 

ground-state  electronic  structure  of  the  compounds.  However, 

comparison  of  the  band  structures  of  Ag  determined  from  photoemission 

and  first-principles  calculations  shows  an  average  deviation  of 

approximately  0.2  eV  in  the  absolute  locations  of  the  valence 

22 

bands,  and  thus  the  photoemission  results  should  represent 
ground-state  bands  to  within  the  uncertainty  of  the  MBBSIS  procedure. 

A  common  method  of  determining  the  quality  of  band  structure 


calculations  is  to  compare  the  ISOS  to  spectral  features  in  valence 
band  XPS  spectra. ^3-25  jn  or^er  to  sake  a  direct  comparison  with  XPS 
spectra,  TDOS  curves  were  produced  with  a  1 .0  eV  Gaussian  broadening 
to  account  for  lifetime  and  instraental  resolution  effects.  The 
overall  agreement  with  the  corresponding  XPS  spectra  vas  good,  as 
shown  in  Fig.  5.  The  detailed  differences  between  the  XPS  spectra 
and  the  TDOS  curves  calculated  from  bands  fitted  to  LAPW  results  may 
be  caused  by  photoemssion  cross  section  effects  or  by  inaccuracies  in 
the  LAPW  bands.  The  TDOS  calculated  using  the  empirically  adjusted 
bands  for  TiN  shows  a  significantly  improved  agreement  with  the  XPS 
spectrom^*  (see  Fig.  5).  For  TiC,  the  agreement  between  the 
empirical  TDOS  and  the  XPS  spectra  was  also  improved  with  respect  to 
the  bands  fitted  to  the  LAPW  results,  but  there  were  still  some 
pronounced  differences.  This  may  have  been  because  fewer 
experimental  ARPES  data1^  for  TiC  were  available  to  use  in  fitting 
the  1BBSIS. 

A  decomposition  of  the  Ti  d-like  TDOS  into  t„  and  e  manifolds 

Zg  g 

provides  interesting  information  about  the  nature  of  the  occupied 
states.  From  the  partial  DOS  calculated  by  the  M3BSIS,  the  member  of 
occupied  Ti  tj^  and  e^  states  is  about  equal  for  TiC,  whereas  there 
are  more  occupied  t-  than  e.  states  in  both  TiN  and  TiO.  Also,  the 
partial  DOS  analysis  showed  that  the  d-like  contribution  to  the 
nonaetal  2p-band  decreased  in  the  series  TiC,  TiN  and  TiO.  This 
indicates  that  the  intermixing  between  noimetal  p  and  metal  d 
orbitals,  or  the  degree  of  covalency,  decreased  going  from  TiC  to 


TiO.  This  agrees  with  the  qualitative  pictures  obtained  f ran  the 


electron  densities.  Blaha,  et  al.1^  compared  the  electron  densities 
of  TiC  and  TiN  calculated  f ran  LAFW  wave  functions  with  the  X-ray 
diffraction  experiments  of  Dunand,  et  al.14  In  this  comparison,  the 
nonspherical  charge  density  around  the  titanimn  atom  determined  from 
the  calculation  and  the  diffraction  data  agreed  very  well. 

Studies  of  chemical  bonding  in  the  Ti  rocksalt  compounds  based 
on  first-principles  calculations  revealed  a  combination  of  ionic, 
covalent  and  metallic  contributions.  Electron  density  plots  may  be 
used  to  show  pictures  of  the  orbital  overlap  or  p-d  interactions,  in 
order  to  visualize  the  chemical  bonding  of  these  compounds.  The 
total  valence  electron  densities  of  TiC,  TiN  and  TiO  calculated  by 
the  1BBSIS  fitted  to  the  first-principles  bands14-16  agree 
qualitatively  with  the  first-principles  results  of  Ref.  9,  even  with 
all  the  approximations  discussed  in  Sec.  II.  The  MBBSIS  electron 
density  plots  for  TiC,  TiN,  and  TiO  exhibit  all  the  same  features  as 
those  shown  in  Fig.  3  of  Ref.  9.  The  electron  densities  around  the 
non-metal  atoms  are  spherically  symmetric  and  show  increased 
localization  going  from  TiC  to  TiO  (Fig.  6).  The  electron  densities 
around  the  titanim  atom  strongly  deviate  from  spherical  symmetry, 
and  the  deviation  changes  from  e^  symmetry  for  TiC  to  t2g  symmetry 
for  TiO. 


The  e^  symmetry  for  the  Ti  electron  density  in  TiC  agrees  with 

9  10 

that  of  Blaha,  et  al.,  but  differs  from  that  of  Trebin  and  Bross. 
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The  latter  charge  density  was  calculated  using  the  Modified  Anptented 

Plane  Vave  (MAW)  method,  and  indicated  a  t^  symmetry  of  the 

electron  charge  density  at  Ti.  This  discrepancy  was  discussed  by 

Blaha,  et  al.,  vho  used  29  k-points  chosen  over  a  uniform  grid  in  the 

irreducible  1/48  of  the  BZ  to  calculate  the  charge  density  in 

Ref.  26.  The  present  study  confirms  that  the  difference  in  the 

symmetry  of  the  charge  density  is  caused  by  the  identity  of  the 

eigenfunctions  used  in  k-simmation  for  calculating  the  total  valence 

electron  densities.  Then  the  total  valence  electron  densities  were 

calculated  with  the  WBSIS  vave  functions  using  the  same  k-points  as 

Trebin  and  Bross,10  the  resulting  plot  shoved  t^  symmetry  at  the  Ti 

atom.  Hovever,  the  total  valence  electron  densities  calculated  using 

1 9 

2,  10  and  60  special  k— points  all  had  e^  symmetry  about  the  Ti 
atom,  and  they  vere  nearly  identical  to  one  another.  Thus,  choosing 
the  correct  k-points  for  calculating  properties  such  as  valence 
electron  densities  is  extremely  important,  but  vhen  using  the  special 
points  of  Chadi  and  Cohen,  only  a  moall  number  of  k-points  are 
required.  The  fact  that  the  MBBSIS  vave  functions  successfully 


reproduced 

the  ch 

am ical 

trends  seen  for 

the 

symmetry 

of  the 

electronic 

charge 

about 

the 

Ti  atom  for 

these 

compounds 

indicate  s 

that  they 

can  be 

useful 

in 

qual itative 

examinations 

of  charge 

densities,  especially  if  first  principles  vave-funct ions  are  not 
available. 

Then  the  total  charge  densities  of  the  16BSIS  fitted  to  LAW 
results  vere  quantitatively  compared  vith  the  first  principles 


resalt*,  significant  differences  were  apparent.  The  differences  were 
particularly  accnte  at  the  aid-point  of  the  Ti-nometal  bond.  In  the 
case  of  TiC,  the  values  of  the  total  charge  densities  aid-way  between 
Ti  and  C  are  8  and  20.7  e/Q  f ran  the  first  principles  resalts  and  the 
M*BSIS,  respectively,  bat  the  aaxiaan  charge  densities  close  to  each 
atca  differ  by  less  than  ~1 0  %.  Thns,  the  1BBSIS  charge  densities 
display  less  localization  than  the  first  principles  resnlts,  which  is 
likely  the  result  of  nsing  snch  a  aall  nunber  of  plane  waves  in  the 
basis.  However,  the  relative  aagnitndes  of  the  NBBSIS  total  charge 
density  aazias  near  Ti  and  the  non-aetal  atca  for  the  three  ccaponnds 

Q 

agreed  very  well  with  the  first  principles  resnlts. 

It  is  also  interesting  to  coapare  the  electron  density  plots  of 
the  eapirically  adjusted  bands  with  those  of  the  LAW  fitted  bands 
for  TiC  and  TiN.^7  By  co* paring  the  NBBSIS  charge  densities  with  each 
other,  the  consequence  of  correcting  the  first  principles  band 
structure  to  agree  with  experiment  should  be  evident.  In  TiC,  the 
eapirically  adjusted  NBBSIS  charge  density  is  higher  at  the  Ti  atom 
and  lower  at  the  C  atca  than  the  NBBSIS  density  directly  fitted  to 
the  LAW  bands.  This  suggests  that  the  experiaental  resnlts  indicate 
aore  covalent  bonding  character  in  TiC  than  predicted  by  the  LAW 
calculations.  In  TiN,  contrary  to  the  TiC  result,  the  adjusted 
NBBSIS  shows  a  lower  electron  density  at  the  Ti  atca  and  higher 
electron  density  at  the  N  at  ton.  Therefore.  TiN  is  saaewhat  aore 
ionic  than  the  LAW  calculation11  indicates.  This  is  coapatible  with 
the  fact  that  the  N-2p  energy  band  deterained  frca  ARPES  data17 
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resides  about  1  eV  farther  below  the  Feral  energy  then  that  of  the 
LAIf  energy  band  and  the  experimentally  det trained  band  dispersions 
are  analler  than  the  calculated  ones. 


An  X-ray  diffraction  experiment  yields  the  total  electron 
density,  bnt  an  energy  band  calculation  can  provide  a  partitioning  of 
electron  density  into  contributions  f ran  individual  bands.  The 
electron  density  plots  of  each  band  for  selected  k-points  can  be  used 
to  visualize  various  types  of  bonding  character,  which  leads  to  an 
improved  understanding  of  the  bonding  properties  of  the  compounds. 


By  comparing  MBBSIS  charge  densities  of  different  compounds, 
qualitatively  correct  chemical  trends  should  be  evident.  The 
relationship  between  banding  and  bonding  may  be  illustrated  by 
examining  states  with  k-points  along  the  A  direction.  There  are  two 
pairs  of  bands  with  the  same  symmetry  Aj  (from  nounetal  2p  and  Ti 
bands)  and  A^  (from  nometal  2p  and  Ti  t2g  bands),  as  well  as  a 
single  A^,  band  (from  Ti  t2g)  and  a  Aj  band  (from  Ti  e^) .  The 
interaction  between  the  states  of  the  same  symmetry  mixes  the  A^ 
states  to  form  p-e^  hybrids  (pda  bonding)  and  the  A^  states  to  form 
P~t2g  hybrids  (pdn  bonding).  In  Fig.  7,  the  electron  density  plots 
for  the  first  band  with  Aj  symmetry  for  the  three  Ti  compounds 
illustrate  the  covalent  bonding  (pdn)  between  nonmetal  2p  and  Ti  t2g 
orbitals.  The  electron  densities  around  the  Ti  atoms  are  reduced 
going  from  TiC  to  TiO,  while  those  around  the  nonmetal  atoms 
increase.  This  indicates  that  TiC  exhibits  the  highest  degree  of 
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covalent  bonding  and  that  ionic  character  increases  as  the  series 
progresses  front  TiC  to  TiO. 

The  electron  density  plots  at  the  Femi  energy  are  illustrated 
in  Fig.  8.  Three  bands  cross  the  Femi  level  in  TiC  and  one  band  in 
TiN  and  TiO  along  the  A  line.  In  TiC,  the  electron  density  plot  of 
the  first  Aj  band  at  the  Femi  energy  reveals  a  slight  cr-type 
covalent  bonding  character  between  the  C  and  Ti  atoms.  In  TiN  and 
TiO,  the  electron  density  plots  of  the  second  band  with  A^  symmetry 
at  the  Femi  level  show  predominantly  metal-metal  interactions 
(Ti  ddQ  sod  4s  in  character).  These  plots  are  important  since  they 

involve  the  electronic  states  that  contribute  most  to  the 

conductivity  of  these  compounds.  Thus,  these  charge  densities 

indicate  that  for  the  perfect  materials  (i.e.  those  for  which 

defects  do  not  dominate  the  charge  transport  properties)  the  Ti  atoms 
would  be  primarily  responsible  for  electron  transport  in  the  Ti 
rocksalt  compounds,  with  the  nometals  acting  essentially  as 
spectators  that  hold  the  crystal  together  by  ionic  forces.  By 
examining  the  charge  density  at  the  Femi  level,  the  importance  of 
Ti-Ti  interactions  becomes  much  more  evident,  since  the  total  valence 
charge  density  is  dominated  by  states  localized  near  the  Ti  and 
nometal  atoms. 

The  disposable  parameters  determined  from  the  fitting  procedure 
for  the  rocksalt  compounds  are  listed  in  Table  I,  along  with  the 
standard  deviation  for  each  parameter.  The  A  parameters,  which  are 
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essentially  three  center  integrals  involving  the  Ti  3d  orbitals,  seen 
to  vary  without  trends  f ran  TiC  to  TiO.  Most  of  the  A  parameters  are 
essentially  zero;  their  absolute  values  are  less  than  their  standard 
deviations.  This  is  expected,  since  the  Ti  atoms  are  separated  by 
fairly  large  distances  in  the  rocksalt  compounds,  and  thus  the  direct 
d-d  interactions  should  be  relatively  nail.  The  much  larger  values 
of  the  Eq  and  A  parameters  in  comparison  with  those  of  the  fee  metals 
might  be  the  result  of  compensating  for  the  large  hybridization 
parameters  Bt  and  Be,  which  are  affected  by  the  tj^  and  e^  band 
positions. 

Unfortunately,  the  parameters  in  the  NBBSIS  do  not  appear  to 
have  a  strong  physical  basis.  Depending  on  the  starting  guesses  in 
the  fitting  procedure,  band  structures  that  are  nearly  identical  can 
arise  from  parameter  sets  with  extremely  different  values.  The  wave 
functions  that  result  from  such  different  parameters  are  also  very 
similar,  and  produce  essentially  identical  charge  densities.  Thus, 
it  appears  as  though  the  NBBSIS  for  rock-salt  type  compounds  is  over 
determined  using  Smith's  parameterization,^  *  and  that  schemes  with 
different  parameters  are  worth  exploring;  in  particular  schemes  that 
involve  only  3  parameters  for  the  direct  d-d  interactions^  and 
additional  pseudopotential  and  orthogonal ization  tens  as  well  as  an 
increased  nimtber  of  plane  waves  should  substantially  improve  the 
NBBSIS  fit  and  make  the  parameters  more  physically  meaningful. 
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IV.  Conclusion 


The  MBBSIS  for  fee  d-band  metals  was  extended  to  inelnde  the 
rocksalt  strnctnre  con  pounds.  The  band  structures  for  TiC,  TiN  and 
TiO  fitted  by  the  present  scheme  were  in  good  agreement  with  the 
first-principles  calculations  up  to  the  tenth  band,  with  an  average 
deviation  of  ~0 .2  eV.  Sem i- empirical  band  structures  for  TiC  and  TiN 
were  also  obtained  by  adjusting  the  interpolated  energy  bands  to 
agree  with  ARPES  data.  This  procedure  assumes  that  the  photoemission 
data  provide  a  good  picture  of  the  ground-state  band  structure.  The 
TDOS  curves  generated  using  the  empirically  adjusted  bands  of  TiC  and 
TiN  demonstrated  better  agreement  with  IPS  valence  band  spectra  than 
the  MBBSIS  bands  fitted  to  the  first-principles  results.  Even  though 
the  approximations  used  for  the  Ti  d-orbital  contribution  to  the 
1BBSIS  wave  functions  were  extreme,  the  valence  electron  density 
plots  calculated  by  using  the  MBBSIS  bands  agreed  qualitatively  very 
well  with  the  first  principles  results  of  Blaha  and  Schwartz. ^  The 
MBBSIS  charge  densities  showed  much  better  agreement  with  those  first 
principles  results  than  did  either  the  EPM  results  of  Ref.  8  or  the 
first  principles  calculation  of  Ref.  10. 

The  1GBSIS  charge  densities  were  used  to  illustrate  several 
points.  The  valence  electron  density  plots  of  the  individual  bands 
showed  an  increasingly  ionic  bonding  contribution  in  the  series  TiC, 
TiN  and  TiO.  The  metal-metal  bonding  is  similar  in  all  these 
compounds,  and  is  most  pronounced  for  energy  bands  at  the  Fenoi 


level.  Cob  pari  ton  of  charge  densities  calculated  from  the  MBBSIS 
fitted  to  first  principles  results  and  the  eapirically  adjusted  bauds 
shoved  that  TiC  involves  sore  covalent  bonding  than  predicted  by  the 
first  principles  calculations,  whereas  TIN  is  sore  ionic. 
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FIGURE  CAPTIONS 


Fig.  1. 

Semi-empirical  energy 

bands 

that 

have 

been 

adjusted 

to 

improve 

agreement  with  the 

ARPES 

data 

of 

Ref. 

17  and 

the 

corresponding  TDOS  cnrve.  The 

energy 

scale 

is 

referenced  to 

the 

Fermi  energy  and  the  TDOS  scale  is 

in  electrons  of  one  spin  per 

eV 

per  unit 

cell  volme. 

Fig.  2. 

Semi-empirical  energy 

bands 

that 

have 

been 

adjusted 

to 

improve 

agreement  with  the 

ARPES 

data 

of 

Ref. 

15  and 

the 

corresponding  TDOS  curve  with  the  sene  remarks  as  Fig.  1. 

Fig.  3.  Interpolated  energy  bands  that  have  been  fitted  to  the 
first-principles  calculation  of  Ref.  16  and  the  corresponding  TDOS 
curve  with  the  same  remarks  as  Fig.  1. 

Fig.  4.  The  comparison  of  the  semi-empirical  and  first-principles 
energy  bands  with  the  ARPES  data  of  Ref.  IS  along  the  A  line  of  the 
BZ  for  TiN.  The  solid  lines  are  the  semi-empirical  energy  bands 
calculated  by  the  1BBSIS,  and  the  dashed  lines  are  the 
first-principles  energy  bands.  The  circles  and  boxes  represent  the 
ARPES  data  of  Ref.  15,  with  the  foxmer  corresponding  to  direct 
transitions  and  the  latter  to  perhaps  indirect  transitions  or  surface 
featnres. 

Fig.  5.  Comparison  between  the  TDOS  corves  of  TiC,  TiN  and  TiO  and 


the  corresponding  ZPS  spectre  in  Ref.  23-25  with  i  piecewise  linear 
background  subtracted.  The  dots  represent  the  adjusted  ZPS  data. 
The  TDOS  carves  of  the  semi-empirical  band  structures  are  shown  as 
solid  lines.  The  TDOS  curves  of  the  interpolated  band  structures 
fitted  to  the  first-principles  calculations  are  represented  as  dashed 
lines.  All  curves  are  plotted  versus  binding  energy  referenced  to 
the  Fensi  energy. 

Fig.  6.  Total  valence  electron  density  plots  for  (a)  TiC.  (b)  TiN 
and  (c)  TiO  in  the  (100)  plane  corresponding  to  the  band  structures 
of  Fig.  1-3,  calculated  by  earning  over  valence  bands  for  10  special 
k-points  and  normalized  to  8,  9  and  10  electrons  per  primitive  unit 
cell,  respectively.  The  contour  lines  represent  intervals  of 
electron  density  of  1  e/0. 

Fig.  7.  Electron  density  plots  of  (a)  TiC,  (b)  TiN  and  (c)  TiO  of 

the  first  band  with  Aj  symmmetry  at  k  *  (2n/a) (0.5 ,0 ,0)  in  the  (100) 

plane.  Each  density  is  normalized  to  2  e/0,  and  the  contour  lines 
represent  intervals  of  electron  density  of  0.5  e/0. 

Fig.  8.  Electron  density  plots  of  (a)  TiC,  (b)  TiN  and  (c)  TiO  for 

the  band  with  Aj  symmetry  at  the  k-point  at  which  the  band  crosses 

the  Fermi  energy  in  Fig.  1-3,  respectively.  Each  density  is 
normalized  to  2  e/0,  and  the  contour  lines  represent  intervals  of 
electron  density  of  1  e/Q. 
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